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ABSTRACT 

Both  numerical  and  experimental  analyses  were  performed 
to  investigate  underwater  shock  propagation  and  the  induced 
nonlinear  response  of  cylindrical  shells  with  end  caps.  The 
cylinders  were  subjected  to  shocks  from  explosive  charges  at 
12  inches  (near-field)  and  28  feet  (far-field)  from  the 
cylinder.  An  underwater  shock  test  was  also  performed  with  the 
far-field  explosion.  The  numerical  results  were  compared  with 
the  experimental  data.  Stresses  and  strains  occurring  in  the 
structure  as  well  as  the  pressure  in  the  water  were  studied. 
The  far-field  explosion  caused  the  largest  circumferential 
deformations  close  to  both  end  plates  and  an  accordion 
oscillatory  motion  of  the  cylindrical  shell.  The  near-field 
explosion  caused  severe  plastic  deformation  in  the 
neighborhood  of  the  closest  end  plate  to  the  charge.  The 
stiffeners  had,  as  expected,  a  larger  effect  on  the 
circumferential  stresses  than  on  the  longitudinal  stresses. 
The  measured  and  calculated  strains  agreed  well  qualitatively 
near  the  remote  end  plate  from  the  charge. 
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I .  INTRODUCTION 

Because  of  the  U.S.  Navy's  high  interest  in  the  underwater 
shock  hardening  effects  on  surface  ships  and  submarines,  this 
research  intends  to  provide  more  insight  into  the  response  of 
a  submerged  vessel  subjected  to  end-on  underwater  shock. 

A  continuous  research  has  been  taken  at  the  Naval 
Postgraduate  School  to  provide  more  insight  into  the 
deformation  and  catastrophic  failure  of  surface  and  subsurface 
hulls.  With  simple  cylindrical  shells  as  a  starting  point,  the 
study  will  then  be  extended  to  structures  with  more  complex 
material  and  geometric  properties  as  the  methods  and 
predictions  improve. 

Some  of  the  previous  studies  in  this  subject  are  listed  in 
references  1-3.  These  studies  have  served  as  building  blocks 
for  the  current  research  into  dynamic  response  of  cylindrical 
shells  to  underwater  shock.  The  objective  of  this  study  is  to 
provide  insight  into  the  end-on  shock  dynamic  response  of 
cylindrical  shells  by  using  numerical  and  experimental 
techniques . 

An  unstiffened  cylinder  subjected  to  a  far-field  explosion 
was  investigated  using  both  numerical  and  experimental 
techniques.  In  addition,  both  unstiffened  and  ring-stiffened 
cylinders  subjected  to  a  near-field  explosion  were  studied 


numerically.  Both  the  finite  element  and  the  boundary  element 
methods  were  utilized  for  the  numerical  study. 


II.  ANALYSIS  AND  EXPERIMENT 

A.   NUMERICAL  ANALYSIS 

The  numerical  study  of  the  behavior  of  cylindrical  shells 
loaded  by  underwater  explosion  was  carried  out  by  using  the 
finite  element  and  boundary  element  methods. 

For  a  problem  with  a  three  dimensional  domain,  the  finite 
element  method  generates  a  three  dimensional  discretization  of 
the  entire  domain  whereas  the  boundary  element  method 
discretizes  the  surface  boundary  of  the  domain  with  a  two 
dimensional  grid.  The  boundary  element  method  reduces 
significantly  the  number  of  elements  required  to  model  the 
problem  by  using  a  two  dimensional  mesh.  Furthermore,  the 
boundary  eloment  method  surpasses  the  finite  element  method  in 
computing  tractions  because  these  tractions  are  treated  as 
primary  but  not  secondary  unknowns  [Ref.  4].  The  matrix 
generated  by  the  boundary  element  method  is  generally  fully 
populated  while  the  matrix  generated  by  the  finite  element 
method  is  usually  narrowly  banded. 

To  study  the  propagation  of  the  explosive  pressure  wave 
through  the  acoustic  medium  and  its  subsequent  interaction 
with  the  cylindrical  shell,  a  finite  element  analysis  program 
called  VEC/DYNA3D  [Ref.  5]  was  used.  This  program  analyzes  the 
deformations  of  solids  which  usually  are  inelastic  in  nature. 


An  alternative  for  modelling  the  acoustic  medium  is  by  using 
the  boundary  element  method  instead  of  using  the  finite 
element  method.  In  doing  so,  the  number  of  elements  is  greatly 
reduced.  The  boundary  element  analysis  method  program  called 
USA  (Underwater  Shock  Analyzer)  [Ref.  6],  was  used  to  compute 
transient  responses  of  submerged  structures  to  acoustic  shock 
waves . 

In  USA,  the  fluid-structure  interaction  was  handled  using 
the  Doubly  Asymptotic  Approximation  (DAA) .  The  differential 
equation  of  motion  describing  the  structure  response  is  given 
below. 

Msx   +  Csx   +  Ksx   =  f  (1) 

where  M£,  C£,  and  KE  are  the  mass,  damping  and  stiffness 
matrices  respectively  and  x,  x,  x  are  the  displacement, 
velocity  and  acceleration  vectors  of  the  structure 
respectively.  The  excitation  for  a  submerged  structure 
subjected  to  an  acoustic  wave  is  given  below. 

f   =  -GAAP,   +  Ps)    +  fd  (2) 

where  -GA,  (P,  +  Ps)  is  the  force  vector  due  to  the  fluid- 
structure  interaction,  fd  is  the  force  vector  applied  to  the 
dry  structure,  P,  is  the  incident  pressure,  Ps  is  the  scattered 
pressure,  G  is  the  transformation  matrix  relating  fluid  and 
structure  nodal  forces  and  Af  is  the  diagonal  area  matrix 
associated  with  the  fluid  mesh. 


To  further  study  the  fluid-structure  interaction  and 
especially  to  relate  the  scattered  wave  pressure  to  velocity 
over  the  wet  surface,  DAA  was  introduced.  The  DAA  solution 
approaches  the  exact  solutions  for  both  early  time  and  late 
time  responses.  The  early  time  response  is  the  high  frequency 
response  and  the  late  time  response  is  the  low  frequency 
response.  The  DAA  is  given  by  the  following  equation. 

MfPs  +  pcAfP£   =  pcMfus  (3) 

where  Mf  is  the  symmetric  fluid  mass  matrix  for  the  wet- 
surface  fluid  mesh,  us  is  the  vector  of  scattered  fluid 
particle  velocities  normal  to  the  structure's  wet  surface,  p 
is  the  density  of  the  fluid,  c  is  the  sonic  speed  in  the  fluid 
and  Af  is  the  symmetric  fluid  mass  matrix  for  the  wet  surface 
fluid  mesh.  This  matrix  is  created  using  the  boundary  element 
method. 

The  high  frequency  approximation,  i.e.  plane  wave 
approximation  implies  that  \P£\  >>\PS\  where  Ps  is  the  time 
derivative  of  the  acoustic  pressure.  Therefore,  Equation  (3) 
is  converted  into  the  following  equation. 

Ps   =  pcu£  (4) 

The  low  frequency  approximation,  i.e.  virtual  mass 
approximation  implies  that  \PS\  <<  \PS\  .  The  concept  of  virtual 
mass   stems   from  the   need  to   include   a  mass   of   fluid 


surrounding  that  structure  at  the  low  frequency  motion  of  the 
structure.  Hence,  Equation  (3)  is  modified  as  seen  below. 

AfPs   =  Mfus  (5) 

Equation  (3)  denotes  the  first  order  of  the  Doubly 
Asymptotic  Approximations  (DAA1) .  The  second  order  of  the 
Doubly  Asymptotic  Approximations,  DAA2,  was  created  to  improve 
the  solution  of  DAA1  for  intermediate  times  and  to  correct  for 
curvature  of  the  surface  of  the  structure  [Ref.  7]. 

The  finite  element  model  of  both  structure  and  fluid  has 
the  advantage  of  presenting  the  pressure  wave  propagation  in 
the  fluid  and  the  fluid-structure  interaction.  On  the  other 
hand,  it  also  requires  a  large  number  of  elements.  Using  the 
finite  element  model  for  the  structure  and  the  boundary 
element  model  for  the  fluid  reduces  the  number  of  degrees  of 
freedom  in  the  system  because  the  fluid  domain  is  usually  much 
larger  than  the  structural  dimension.  However,  it  cannot 
present  the  propagation  of  the  pressure  wave  across  the  fluid. 

The  post-processing  of  the  VEC/DYNA3D  and  USA  programs  was 
done  using  LS-TAURUS  [Ref.  8].  This  post-processor  generates 
pressure,  strain  and  stress  contours  superposed  on  the  mesh  or 
also  generates  element  and  node  time  history  responses. 

B.   NUMERICAL  MODELLING 

To  study  the  deformation  of  a  cylindrical  shell  subjected 
to  an  end-on  shock,  two  numerical  models  were  created.  The 


first  model  was  designed  to  study  the  propagation  of  the 
pressure  wave  from  the  explosive  to  the  cylinder  and  its  early 
interaction  with  the  cylindrical  shell.  The  second  model  was 
designed  to  study  stresses  and  strains  of  the  cylindrical 
shell  without  analyzing  the  fluid  or  explosive  around  it. 

To  generate  the  first  model,  the  pre-processor ,  INGRID 
[Ref.   9],   generated  the   finite  element  meshes   for  the 
explosive,  fluid  and  cylindrical  shell  and  VEC/DYNA3D  computed 
their  dynamic  response. 

The  computational  effort  was  minimized  by  creating  a 
quarter  model,  possible  only  because  of  the  symmetric  geometry 
of  the  charge  and  cylindrical  shell  locations.  Appropriate 
symmetric  boundary  conditions  were  applied  to  the  problem. 

The  explosive  was  modeled  with  a  fine  mesh  of  416  elements 
to  avoid  non-spherical  propagation  and  severe  distortion  in 
the  finite  element  mesh  as  a  result  of  the  expansion  of  the 
explosive  in  the  fluid. 

To  model  the  explosion,  the  Jones-Wilkins-Lee  (JWL) 
equation  of  state  was  invoked  to  describe  the  pressure-volume- 
energy  behavior  of  high  explosives.  [Ref.  10] 
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where  A,  B,  and  C  are  linear  coefficients  given  in  units  of 
Mbar.  Rl,  R2  and  w  are  nonlinear  coefficients,  V  is  the  volume 
of  detonation  products  divided  by  the  volume  of  undetonated 


high  explosive,  P  is  the  pressure,  and  E  is  the  detonation 
energy  per  unit  volume  in  (Mbar  cm3) /cm3. 

The  Gruneisen  equation  of  state  was  used  to  define  the 
pressure  for  compressed  materials,  in  this  case  water.  This 
equation  of  state  is  provided  in  [Ref.  10]. 

The  model  including  the  charge,  water  and  cylinder  had 
11808  elements  which  consequently  required  large  amounts  of 
storage.  Finite  element  meshes  are  given  in  Figure  2.1. 
Therefore,  the  charge  was  constrained  to  be  close  to  the 
target  cylinder  to  minimize  the  elements  between  explosive  and 
cylinder  and  still  be  able  to  store  the  shock  wave  propagation 
and  stress  wave  effect  information.  The  cylindrical  shell  of 
this  finite  element  model  was  modified  to  add  ring  stiffeners. 
One  model  had  one  ring  stiffener  located  halfway  between  the 
cylinder  end  plates.  The  second  model  had  two  ring  stiffeners 
equidistant  between  the  two  end  plates.  A  view  of  these  two 
models  is  shown  in  Figure  2.2.  The  additional  stiffeners  only 
added  an  additional  16  elements  per  stiffener  to  the  total. 
The  ring  stiffeners  were  added  to  study  their  effect  in  stress 
wave  propagation  and  reinforcement  properties. 

The  other  model  used  in  this  research  was  designed  to 
study  the  stresses  and  strains  of  the  cylinder  for  comparison 
with  experimental  results  of  an  actual  underwater  shock  test. 
The  model  dynamic  response  was  analyzed  using  both  VEC/DYNA3D 
and  USA.  The  finite  element  method  was  used  for  the 
cylindrical  shell  and  the  boundary  element  method  was  used  for 
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(d)  Quarter  Fluid  Mesh 
Figure  2.1.  Finite  Element  Meshes  for  Near-Field  Explosion 


(a)  One-Stif fener  Model 


(b)  Two-Stif fener  Model 
Figure  2.2.  Stiffened    Cylinders    with    End    Plates 
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the  fluid-structure  interaction.  This  model  consisted  of  a 
cylindrical  shell  capped  at  both  ends  by  two  end  plates.  The 
cylindrical  shell  had  512  elements  and  the  endplates  had  48 
elements  each.  The  full  model  had  608  total  elements.  A 
diagram  of  this  model  is  presented  in  Figure  2.3.  The 
different  size  of  the  cylindrical  shell  elements  was  generated 
to  get  a  better  solution  in  the  strain  gage  locations  of  the 
experiment.  These  locations  will  be  presented  when  the 
experimental  procedure  is  explained. 

A  study  of  grid  independence  and  time  step  instability  was 
performed  to  verify  that  the  choice  of  the  element  sizes  and 
time  steps  chosen  were  not  affecting  the  solution.  To  analyze 
grid  independence,  three  different  discretizations  were  made 
over  a  four  inch  segment  of  the  cylinder.  The  discretization 
was  made  in  the  longitudinal  direction  due  to  the  symmetry  of 
the  problem.  The  coarse  model  had  two  elements  of  two  inches 
each  in  the  axial  direction.  The  finer  model  had  six  elements 
of  0.66  inches  each  in  the  same  direction.  The  finest  mesh 
model  had  eighteen  elements  of  0.22  inches  each.  The  three 
discretizations  are  shown  in  Figure  2.4.  Comparison  of  the 
three  models  revealed  considerable  differences  between  the 
coarse  model  and  the  two  finer  mesh  models.  The  six  and  18 
element  models  had  almost  identical  results.  As  a  result  of 
this  comparison,  the  six  element  model  was  chosen  as  a  good 
estimate  of  the  size  of  the  elements  to  use  for  the  actual 
model . 
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(b)  Cylindrical  shell 

Figure  2.3.  Finite  element  meshes  of  a  cylinder  subject  to  a 
far-field  explosion 
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(a)  Coarse  model  (2  inch  elements) 
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(b)  Finer  model  (0.66  inch  elements 
(6  elements  in  4  inch  span) 
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(c)  Finest  model  (0.22  inch  elements) 
(18  elements  in  4  inch  span) 

Figure  2.4.  Finite  element  discretization  over  4  inch  span  at 
one  end  of  cylinder 
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Another  area  of  concern  was  the  effect  the  time  step  would 
have  on  the  solution.  An  appropriate  method  to  choose  a  stable 
time  step  size  was  the  Courant-Friederichs-Lewy  Criterion.  The 
equation  for  this  criterion  is  given  below. 

At<C/c  (7) 

where  At  is  the  time  step  size,  d      is  the  length  of  the 
element,  and  c  is  the  sonic  speed  in  6061-T6  Aluminum. 

Two  different  time  step  sizes  were  chosen,  2xl0e~6  and 
2xl0e~7  seconds.  These  time  steps  satisfied  the  above  mentioned 
criterion.  The  resulting  solutions  for  the  two  time  steps  were 
almost  identical.  Therefore,  4xl0e~7  seconds  was  chosen  as  the 
time  step  for  this  problem  since  this  value  was  between  the 
two  chosen  for  the  criterion  verification. 

C.   UNDERWATER  EXPLOSION  TEST 

The  underwater  explosion  test  was  performed  at  Dynamic 
Testing  Inc.  (DTI)  facilities  in  Rustburg,  Virgina.  This 
facility  had  a  quarry  that  had  been  filled  with  water  for  use 
in  underwater  shock  tests.  The  water  depth  was  approximately 
130  feet  at  the  location  of  the  test  which  was  deep  enough  to 
allow  the  study  of  the  dynamic  response  of  the  cylinder  prior 
to  the  arrival  of  the  reflected  shock  wave  from  the  bottom. 

The  cylinder  was  placed  12  feet  below  the  water's  surface 
and  held  in  place  by  a  crane  with  pendants  attached  to  both 
front  and  rear  end  plates.  The  rig  attachment  can  be  observed 
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in  Figure  2.5.  The  60  pound  charge  of  HBX1  was  also  placed  12 
feet  below  the  surface  and  aligned  with  the  cylinder  with  a 
span  wire  from  the  charge  float  to  the  crane  rig.  Figure  2.6. 
gives  a  good  view  of  the  arrangement.  The  12  foot  depth  was 
chosen  so  that  the  bubble  generated  by  the  explosion  would 
vent  to  the  surface  prior  to  encountering  the  cylinder.  The 
test  geometry  is  shown  in  Figure  2.7.  The  explosive  charge  was 
activated  by  a  radio  device  and  the  plume  is  pictured  in 
Figure  2.8. 

The  strain  gages  used  for  this  test  were  of  type  CEA-06- 
250UW-350.  These  are  general  purpose  strain  gages  with  an 
optimum  operating  range  of  ±  1500  micro  strain  and  are  used 
for  both  static  and  dynamic  test  measurements.  They  were 
attached  to  the  cylinder  using  a  M  bond  200.  There  were  a 
total  of  seven  strain  gages  placed  at  locations  A,  B  and  C,  as 
seen  in  Figure  2.9,  per  axis  for  a  total  of  14  strain  gages. 
The  cylinder  was  oriented  so  that  the  gages  at  C  would  be 
closest  to  the  explosive  charge.  A  pressure  probe  was  also 
placed  28  feet  from  the  cylinder  to  measure  the  free  field 
pressure . 

The  cylinder  used  for  the  underwater  shock  test  was 
manufactured  from  6061-T6  Aluminum.  This  alloy  is  primarily  an 
Aluminum-Magnesium-Silicon  alloy.  The  T6  denomination 
indicates  it  was  solution  treated  and  artificially  aged  [Ref. 
11].  The  cylinder  consisted  of  a  quarter  inch  thick 
cylindrical  shell  and  two  one-inch  thick  circular  end  plates 
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Figure  2.5.  Cylinder  and  Crane  Rigging 
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Figure  2.6.  Test  Geometry  at  the  Surface 
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Figure  2.7.  Test  Geometry  (Ranges  in  Feet) 


Figure  2.8.  Water  Plume  from  Shot 
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Figure  2.9.  Cylinder  Geometry 
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as  seen  in  Figure  2.9.  The  welding  of  the  two  end  plates  to 
the  cylindrical  shell  was  done  at  the  Naval  Postgraduate 
School  facilities.  There  was  a  concern  in  the  welding  process 
of  the  two  circular  end  plates  to  the  shell  because  of  the 
heat  generated  at  the  weld  and  its  effect  on  the  alloy 
morphology  close  to  the  weld.  However,  the  strain  readings 
were  taken  far  enough  from  the  heat  affected  zone  to  have  any 
measurable  effect  on  the  readings.  The  welding  was  done  using 
tungsten  inert  gas  (TIG)  .  This  procedure  has  been  recommended 
by  most  expert  welders  [Ref .  12] . 

The  cylinder  weighed  60.5  pounds  and  tensile  tests  were 
done  to  verify  that  the  material  properties  were  close  to  the 
nominal  properties  of  6061-T6  Aluminum.  This  test  can  be  seen 
in  the  Appendix.  The  value  of  Young's  Modulus  was  10800  ksi 
and  the  yield  strength  was  approximately  43  ksi. 

Post-shot  observation  showed  no  visible  deformations  on 
the  cylinder  as  seen  in  Figure  2.10.  All  strain  gages  were 
well  fixed  when  uncovering  the  bonding  material.  Only  strain 
gage  CI  in  Figure  2.9.  was  found  wet  when  uncovered  suggesting 
the  possibility  of  water  insertion. 


21 


u 

•8 

a 

■H 
H 
>i 
O 

(H 
O 

0) 

•rl 
> 

•P 

o 

& 
w 

i 
■p 

0) 

o 

04 


CM 

<D 
U 

3 

o> 

•H 
fa 


22 


III.  RESULTS  AND  DISCUSSIONS 

A.   NUMERICAL  RESULTS  OF  FAR-FIELD  EXPLOSION 

A  cylinder  was  subjected  to  an  end-on  explosion  with  a 
standoff  distance.  The  sketch  of  the  configuration  was 
presented  in  previous  Figure  2.6.  The  deformations  of  the 
cylinder  are  shown  in  Figure  3.1  at  different  times  as  the 
shock  wave  propagates  along  the  axial  direction  of  the 
cylinder.  Because  of  the  symmetric  loading,  the  deformation 
was  axially  symmetric.  In  order  to  visualize  the  deformations 
more  clearly,  the  actual  deformations  were  magnified  by  the 
scale  factor  shown  in  Figure  3.1. 

As  the  shock  pressure  wave  hit  one  end  plate,  which  is  at 
the  nearest  location  from  the  charge,  the  compressive  stress 
wave  propagated  from  the  end  plate  to  the  other  end  plate 
through  the  aluminum  alloy  cylinder  with  a  faster  speed  of 
sound  than  the  shock  pressure  wave  propagating  through  the 
surrounding  water  medium.  The  nearer  end  plate  will  be  called 
the  near  plate  while  the  other  end  plate  will  be  called  the 
remote  end  plate  in  the  following  discussion.  The  shock  load 
caused  initially  a  localized  circumferential  deformation  close 
to  the  near  end  plate  as  shown  in  Figure  3.1c.  This  localized 
deformation  remained  as  a  permanent  deformation.  As  the  shock 
pressure  wave  propagated  toward  the  remote  end  plate,   it 
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induced  a  radial  contraction  of  the  cylinder  wall  and 
subsequently  a  recovery  of  it.  After  the  shock  pressure  wave 
passed  the  remote  end  plate,  there  was  also  a  localized 
circumferential  deformation  near  the  remote  end  plate  as  shown 
in  Figure  3.1h.  The  overall  steady  state  deformation  of  the 
cylinder  was  nearly  symmetric  about  the  center  plane  which  is 
located  at  an  equal  distance  between  the  two  end  plates.  At 
this  stage,  the  cylinder  had  the  deformation  near  the  two  end 
plates,  and  the  rest  of  the  cylinder  had  very  little  radial 
deformation.  The  axial  deformation  of  the  cylinder  had  an 
accordion  mode  as  expected.  Figure  3.2  is  the  plot  of  axial 
velocities  at  the  centers  of  the  two  end  plates.  The  two 
velocities  had  a  phase  difference  of  180  degrees.  The  phase 
difference  indicated  the  accordion  mode. 

The  initiaJ  localized  circumferential  deformation  near  the 
end  plates  may  be  explained  as  follows.  The  axisymmetric 
deformation  of  a  shell  has  the  following  governing  equation. 

D^w+F^w+Eh_^=P  (8) 

dx<  dx2  a- 

where  w  is  the  radial  deflection,  x  is  the  axial  direction,  D 
is  the  flexural  rigidity  of  the  shell,  E  is  the  elastic 
modulus,  h  is  the  wall  thickness  of  the  shell,  a  is  the  shell 
radius,  and  F  and  p  are  the  axial  load  and  the  radial  load, 
respectively.  Therefore,  the  effect  of  the  inertia  term  on  the 
deformation  was  negligible.  The  axial  load  applied  on  the 
cylinder  due  to  the  shock  pressure  was  lower  than  the  critical 
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buckling  load  of  that  cylinder  from  the  linear  buckling 
theory.  In  this  case,  the  deformation  of  the  cylinder  was  of 
a  decaying  exponential  form  and  it  damped  out  rapidly  if  the 
cylinder  was  not  short  [Ref .  13] .  The  very  stiff  end  plates 
suppressed  the  deformation  very  near  the  end  plates.  This 
resulted  in  a  shift  of  the  locations  of  the  circumferential 
deformation  somewhat  away  from  the  end  plates.  The  inertial 
force  was  not  included  in  the  equation.  When  the  density  of 
the  cylinder  was  varied  from  its  nominal  value  to  a  value  of 
one  order  of  less  magnitude,  the  same  kind  of  initial 
localized  deformation  was  observed.  However,  the  low  density 
caused  less  permanent  localized  deformation. 

Two  normal  strain  components,  i.e.,  hoop  (or 
circumferential)  and  axial  (or  longitudinal)  strains,  were 
computed  at  some  selected  locations.  The  selected  points  were 
located  at  the  locally  deformed  zones  as  well  as  at  the  center 
of  the  cylinder.  The  locations  where  the  strains  were  computed 
are  illustrated  in  Figure  3.3.  Strain  gages  were  also  attached 
to  the  same  locations  on  the  cylinder  in  the  experimental 
study  for  the  comparison  with  the  numerical  study.  All  the 
strains  were  computed  and  measured  at  the  outside  surface  of 
the  cylinder.  The  comparison  between  the  numerical  and 
experimental  results  will  be  provided  in  the  next  section.  As 
shown  in  the  strain  plots  of  Figure  3.4,  the  circumferential 
strains  at  locations  A  and  C  of  Figure  3.3  were  very  close  to 
each  other.  The  strains  had  initially  large  compressive  peaks 
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and  returned  to  the  steady  state  values  which  were  about  half 
the  initial  peak  strain  values.  Location  A  had  a  small 
circumferential  tensile  strain  before  arrival  of  the  shock 
pressure  wave.  The  small  tensile  strain  was  caused  by  the 
longitudinal  compressive  stress  wave  which  arrived  at  the 
location  through  the  cylinder.  The  longitudinal  compressive 
strain  caused  by  the  longitudinal  stress  wave  resulted  in  the 
circumferential  tensile  strain  because  of  Poisson's  effect. 
After  the  shock  pressure  wave  arrived  at  the  location,  the 
circumferential  strain  became  compressive.  The  circumferential 
strain  at  the  middle  of  the  cylinder,  i.e,  location  B  in 
Figure  3.3,  was  much  less  that  those  at  locations  A  and  C  as 
shown  in  Figure  3.5  because  the  localized  deformations  at 
locations  A  and  C  induced  the  larger  circumferential  strains. 
The  initial  peak  of  longitudinal  strain  was  compressive  at 
all  locations  due  to  the  propagation  of  the  longitudinal 
compressive  stress  wave  after  the  shock  pressure  hit  the  near 
end  plate.  The  early  arrival  of  the  shock  pressure  wave  to 
location  C  in  Figure  3.3,  caused  the  immediate  jump  of  the 
longitudinal  strain  into  tension  after  the  initial 
compression.  On  the  other  hand,  the  subsequent  compressive 
stress  wave  resulted  in  more  compression  at  locations  A  and  B 
than  at  location  C  until  the  shock  pressure  wave  arrived  at 
the  locations.  At  later  times  the  longitudinal  strain  were 
quite  different  at  locations  A  and  C.  The  longitudinal  strain 
at   location  A  stayed  in  tension   just  after  the  initial 
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Figure    3.5. 


Circumferential  Strain  20.75  Inches  from 
Closest   End  Plate  to  the  Explosive  Charge 
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compressive  peak  while  the  longitudinal  strain  at  location  C 
alternated  between  tension  and  compression.  However,  the 
longitudinal  strain  was  quite  sensitive  with  the  locations 
near  A  and  C  because  of  the  localized  deformation. 
Longitudinal  strains  at  various  locations  near  position  C  were 
plotted  in  Figure  3.6.  The  approximate  frequency  of  the 
longitudinal  strains  after  the  initial  compression  was  about 
750  Hz  at  all  locations. 

The  present  study  showed  that  the  end  plates  played  an 
important  role  in  the  deformation  of  the  cylinder.  Some 
parametric  study  was  performed  to  find  the  effect  of  the  end 
plates  on  the  dynamic  response  of  the  cylindrical  shell.  The 
first  study  was  to  find  out  the  effect  of  the  inertial  force 
of  the  end  plates.  Therefore,  the  density  of  the  end  plates 
was  reduced  tenfold  with  the  stiffness  of  the  plates  remaining 
the  same.  The  light  end  plates  had  a  larger  effect  on  the 
longitudinal  strains  than  on  the  circumferential  strains.  The 
frequency  of  the  longitudinal  strains  was  approximately  840  Hz 
for  the  light  end  plates.  This  frequency  was  higher  than  that 
for  the  heavy  plates  as  expected.  The  longitudinal  strain 
damped  more  quickly  with  the  light  end  plates  than  with  heavy 
end  plates  as  shown  in  Figure  3.7.  The  light  end  plates 
resulted  in  an  approximately  30%  decrease  of  the 
circumferential  strain  at  location  C  while  its  effect  was 
negligible  on  the  circumferential  strain  at  location  A. 
Moreover,  the  circumferential  strain  at  location  B  increased 
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A.  (165)  4.0  inches  from  end  plate 

B.  (169)  4.5  inches  from  end  plate 

C.  (173)  5.0  inches  from  end  plate 
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Figure  3.6. 


Longitudinal   Strains  in  the  Vicinity  of  the 
Closest  End  Plate  to  the  Explosive  Charge 
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about  20  percent  for  the  light  end  plates  as  shown  in  Figure 
3.8. 

The  next  parametric  study  was  performed  to  investigate  the 
effect  of  stiffness  of  the  end  plates  on  the  deformation  of 
the  cylinder.  The  elastic  modulus  of  the  plate  was  reduced 
tenfold  without  change  of  the  rest  of  the  material  data.  The 
less  stiff  end  plates  made  an  effect  on  not  only  longitudinal 
strain  but  also  circumferential  strain  as  shown  in  Figures  3.9 
and  3.10.  The  less  stiff  support  of  the  cylinder  at  the  end 
plates  caused  more  fluctuation  in  the  circumferential  strain 
as  shown  in  Figure  3.9.  The  wave  pattern  in  the  longitudinal 
strain-time  history  plot  was  also  severely  altered  by  the  less 
stiff  end  plates  as  shown  in  Figure  3.10.  The  longitudinal 
velocities  of  the  centers  of  the  less  stiff  end  plates  were 
plotted  in  Figure  3.11  and  showed  a  very  different  pattern 
than  the  accordion  oscillation  of  the  stiff  plates  in  Figure 
3.2.  The  density  change  in  the  previous  parametric  study  did 
no  show  a  different  velocity  pattern  of  the  end  plates. 

B.   COMPARISON  BETWEEN  NUMERICAL  AND  EXPERIMENTAL  RESULTS 

An  underwater  explosion  test  was  performed  at  Dynamic 
Testing  Inc.  facilities  in  Rustburg,  Virginia  as  described  in 
Chapter  II.  Unfortunately,  a  pressure  gage,  placed  to  measure 
a  free-field  pressure  at  the  same  stand-off  distance  as  the 
cylinder,  failed.  Therefore,  no  information  was  available  for 
the  shock  pressure  due  to  the  charge.  The  previous  numerical 
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Figure  3.7.  Longitudinal  Strains  37.5  Inches  from  the  Closest 
End  Plate,  of  Different  Densities,  to  the 
Explosive  Charge 
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(a)  nominal  density  end  plate 
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(b)  lower  density  end  plate  (order  of  ten) 


Figure  3.8. 


Circumferential  Strains  20.75  Inches  from  the 
Closest  End  Plate,  of  Different  Densities,  to  the 
Explosive  Charge 
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(a)  nominal  stiffness  end  plate 
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(b)  lower  stiffness  end  plate  (order  of  ten) 


Figure  3.9 


Circumferential  Strains  37.5  Inches  from  the 
Closest  End  Plate,  of  Different  Stiffness,  to  the 
Explosive  Charge 
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(a)  nominal  stiffness  end  plate 
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Figure  3.10. 


Longitudinal  Strains  37 . 5  Inches  from  the 
Closest  End  Plate,  of  Different  Stiffness,  to 
the  Explosive  Charge 
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Figure  3.11 


Nodal  Velocities  of  Both  End  Plates  of 
Stiffness  Ten  Times  Lower  Than  the  Nominal 
Value  for  6061-T6  Aluminum 
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simulation  was  carried  out  using  the  empirical  equation  for 
pressure  with  the  nominal  weight  of  the  charge.  As  a  result, 
the  following  comparison  was  made  in  a  qualitative  sense.  In 
addition,  several  strain  gages  failed  before  and  during  the 
test.  All  the  strain  gages  to  measure  circumferential  strains 
failed.  Only  the  axial  strains  were  compared  here. 

The  speed  of  sound  in  the  water  was  computed  based  on  the 
arrival  time  of  the  pressure  wave  to  the  strain  gage  from  the 
explosion  and  the  stand-off  distance.  The  computed  speed  of 
sound  was  close  to  a  nominal  value  of  5000  ft/sec.  Then 
locations  of  the  strain  gages  were  shown  in  Figure  3.3.  The 
longitudinal  strain  at  gage  location  A  was  compared  in  Figure 
3.12.  All  experimental  data  was  filtered  out  at  2000  Hz  low 
pass.  Therefore,  there  was  no  peak  strain  with  a  higher 
frequency  than  2000  Hz  in  the  experimental  data.  Both 
numerical  and  experimental  solutions  indicated  the  initial 
compression  and  tension  at  later  times.  They  agreed  well 
qualitatively  even  if  there  was  a  mismatch  in  magnitude.  The 
longitudinal  strain  was  quite  sensitive  to  the  location.  The 
longitudinal  strain  computed  at  just  one  inch  away  from  the 
gage  location  was  compared  with  the  experimental  measurement 
in  Figure  3.13.  Two  strain  gages,  both  of  which  were  located 
at  location  A  but  were  separated  in  the  circumferential 
direction  with  an  angle  of  180  degrees,  gave  a  similar 
longitudinal  measurement  as  shown  in  Figure  3.12.  This 
measurement   indicated   the   axisymmetric   nature   of   the 
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(a)  AIL  Strain  Gage  Location  (See  Figure  2.9) 


Figure  3.12 


Measured  and  Computed  Longitudinal  Strains  at 
37.5  Inches  from  Closest  End  Plate  to 
Explosive 
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(b)  A3L  Strain  Gage  Location  (See  Figure  2.9) 


Figure  3.12 


(cont.)  Measured  and  Computed  Longitudinal 
Strains  at  37.5  Inches  from  Closest  End  Plate 
to  Explosive 
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deformation.  However,  the  two  strain  gages,  which  were  placed 
at  location  C  and  were  separated  180  degrees  in  the 
circumference,  provided  quite  different  longitudinal  strains. 
This  discrepancy  in  readings  could  have  been  caused  due  to 
observed  water  insertion  in  the  strain  gage  at  location  CI. 

Accordingly,  the  numerical  result  was  not  in  good 
agreement  with  the  experimental  data  as  shown  in  Figure  3.14. 
The  sensitivity  of  the  location  were  the  strains  were  computed 
was  also  evident  for  location  C  as  shown  in  Figure  3.15.  The 
experiment  showed  a  compressed  longitudinal  strain  while  the 
numerical  solution  showed  an  oscillation  of  the  strain  from 
compression  to  tension.  A  small  distance  away  from  location  C 
gave  a  different  solution.  However,  the  experimental  data  was 
still  different  from  the  numerical  result.  The  longitudinal 
strain  at  gage  location  B  was  compared  in  Figure  3.16.  The 
result  was  similar  to  that  at  gage  location  C.  It  was  not 
clear  at  this  time  what  the  major  cause  of  the  discrepancy 
was.  However,  this  was  the  first  test  of  a  series  of 
underwater  experiments  to  be  performed.  The  discrepancy  will 
be  investigated  from  the  following  experiments  with  more 
available  experimental  data  and  numerical  simulations. 

C.   NUMERICAL  RESULTS  OF  NEAR-FIELD  EXPLOSION 

The  study  of  the  dynamic  response  of  cylindrical  shells 
subjected  to  a  close-in  explosion  was  performed.  The  stand-off 
distance  was  12  inches  from  one  end  plate.   The  overall 
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Figure  3.13 


Measured  and  computed  longitudinal  strains  at 
38 . 5  inches  from  closest  end  plate  to 
explosive 
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(a)  C1L  Strain  Gage  Location  (See  Figure  2.9) 


Figure  3.14 


Measured  and  Computed  Longitudinal  Strains  at 
4.5  Inches  from  Closest  End  Plate  to 
Explosive 
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(b)  C3L  Strain  Gage  Location  (See  Figure  2.9) 


Figure  3 . 14 


(cont.)  Measured  and  Computed  Longitudinal 
Strains  at  4 . 5  Inches  from  Closest  End  Plate 
to  Explosive 
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Measured  and  Computed  Longitudinal  Strains  at 
3.5  Inches  from  Closest  End  Plate  to 
Explosive 
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B2L  Strain  Gage  Location  (See  Figure  2.9) 


Figure  3.16 


Measured  and  Computed  Longitudinal  Strains  at 
20.75  Inches  from  Closest  End  Plate  to 
Explosive 
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geometry  with  finite  element  meshes  is  shown  in  Figure  2.1.  A 
two  pound  explosive  charge  of  pentolite  was  used  in  the 
numerical  study.  The  dimension  of  the  cylinder  was  the  same  as 
for  the  far-field  explosion  except  for  a  length  of  43  inches 
instead  of  42  inches  long  and  the  addition  of  ring  stiffeners. 
1 .  Pressure  Wave  Study 

The  explosive  charge  was  refined  in  the  numerical 
model  to  avoid  discontinuities  of  the  sphere  as  it  expanded 
with  time.  The  pressure  wave  propagated  spherically  through 
the  finite  element  mesh  as  shown  in  Figure  3.17. 

The  wave  propagation  was  studied  especially  at  the 
first  interaction  between  the  pressure  wave  and  one  end  plate. 
The  cylinder  was  accelerated  when  the  wave  first  interacted 
with  the  end  plate.  At  about  0.3  milliseconds  the  fluid  and 
end  plate  velocities  break  away  for  a  short  period.  The  fluid 
velocity  was  lower  than  that  for  the  shell.  Consequently, 
tension  would  be  induced  in  the  fluid  but  since  fluids  cannot 
experience  tension,  the  water  particles  break  away  creating  a 
vacuum.  This  effect,  known  as  hull  cavitation,  can  be  observed 
in  Figure  3.18. 

2  .  Ring  Stif fener  Study 

The  three  cylinders  studied  were  an  unstiffened 
cylinder,  a  one-ring  stiffener  cylinder,  and  a  two-ring 
stiffened  cylinder.  The  ring  stiffeners  were  located  at  equal 
intervals  between  the  end  plates.  In  other  words,  one  ring 
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Figure  3.17 


Pressure  Wave  Propagation  History  in  Near- 
Field  Explosion  Model 
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Figure  3.18 


Cavitation  at  End  Plate  Closest  to  Near-Field 
Explosion 
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stiffener  was  placed  at  the  center  of  the  cylinder  as  seen 
previously  in  Figure  2.2a  and  the  two  ring  stiffeners  were 
placed  at  one  third  distances  from  the  end  plates  respectively 
as  seen  previously  in  Figure  2.2b.  The  stiffeners  were  one 
inch  high  and  0.25  inches  thick.  The  nominal  material  property 
data  of  a  6061-T6  Aluminum  alloy  was  used  for  this  study. 

The  ring  stiffened  cylinders  were  compared  to  the 
unstiffened  cylinder  to  provide  insight  into  the  effectiveness 
of  the  ring  stiffeners  in  affecting  stress  wave  propagation 
and  deformation  of  cylinders.  The  numerical  model  was  labeled 
in  specific  locations  where  the  dynamic  response  was  of 
interest.  These  locations  are  indicated  in  Figure  3.19  and 
they  include  the  largely  deformed  area  at  the  front  end  of  the 
cylinder,  F,  the  stiffener  locations:  SI,  S2,  and  S3,  and  the 
location  close  to  the  remote  end  plate,  R. 

The  initial  shock  pressure  produced  a  severe  local 
deformation  about  three  inches  from  the  front  end  plate  for 
all  three  cases  of  cylinders  as  shown  in  Figure  3.20.  The 
cylinder  yielded  identically  in  all  cases  because  the 
stiffeners  had  no  effect  on  the  initial  deformations  near  the 
front  end  plate  until  the  pressure  wave  passed  through  the 
stiffeners.  Comparison  of  the  circumferential  strains  at 
location  SI  between  the  unstiffened  and  two-stif fener  models 
was  given  in  Figure  3.21.  In  the  figure  legends  A  and  B 
denoted  shell  elements  located  just  before  and  after  the 
stiffener.  The  stiffener  reduced  the  compressive  peak  strain 
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at  the  location  of  the  stiffener  as  expected.  However,  the 
unstiffened  cylinder  had  a  larger  relief  of  the  compressive 
strain  than  the  two-stiffened  cylinder.  Comparing  the  same  at 
location  S3  indicated  that  the  stiffener  not  only  reduced  the 
compressive  peak  strain,  but  also  relieved  the  strain  more  at 
later  times  as  shown  in  Figure  3.22.  The  circumferential 
strains  at  the  two  stiffeners  were  plotted  in  Figure  3.23.  The 
front  stiffener  had  a  much  larger  strain  than  the  back 
stiffener.  The  large  plastic  strain  at  the  front  stiffener 
restrained  the  recovery  of  the  circumferential  strain  at  the 
adjacent  shell  elements  at  later  times  compared  to  the 
unstiffened  cylinder  as  shown  in  Figure  3.22. 

Comparison  of  the  unstiffened  and  one-stif fener 
cylinders  revealed  that  the  stiffener  caused  a  large  relief  of 
the  circumferential  strain.  The  center  stiffener  also  caused 
a  big  difference  in  the  circumferential  strains  between  the 
shell  elements  just  before  and  after  the  stiffener.  The  strain 
was  greatly  reduced  after  the  stiffener  as  seen  in  Figure 
3.24.  The  longitudinal  strains  were  compared  in  Figures  3.25 
and  3.26  between  two-stiffened  and  unstiffened.  Figure  3.25 
was  the  comparison  at  location  SI  and  Figure  3.26  at  location 
S3.  The  ring  stiffener  at  location  SI  altered  the  longitudinal 
strain  before  and  after  the  stiffener  as  shown  in  Figure  3.25. 
The  stiffener  at  location  SI  induced  a  higher  longitudinal 
strain  at  the  shell  element  located  just  before  the  stiffener 
compared  to  that  at  the  same  location  of  the  unstiffened 
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(a)  Unstiffened  model 


(b)  One-stif fener  model 


(c)  Two-stif fener  model 
time  =  .35999E  +  03  |is 

Figure  3.20.     Local  Deformations  for  Near-Field  Models 
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Figure  3.21. 


Circumferential  Strains  at  Both  Sides  of 
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(b)  One-stif fener  model 


Figure  3.22 


Circumferential  Strains  at  Both  Sides  of 
Stiffener  Farthest  from  Explosive  Charge 
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(b)  Stiffener  28.5  in.  from  end  close  to  explosive 


Figure  3.23 


Circumferential  Strains  at  Both  Stiffeners  of 
Two-Stiff ener  Model 
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cylinder.  In  addition,  the  stiffener  relieved  the  strain  near 
to  zero  at  the  shell  element  located  just  after  the  stiffener 
as  time  elapsed.  The  effect  of  the  stiffener  at  location  S3  on 
the  longitudinal  strain  was  much  smaller  compared  to  that  at 
location  SI.  No  accordion  mode  was  observed  for  all  three 
cases  of  cylinders. 
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(a)  Unstiffened  model 


(b)  One-stif fener  model 


Figure  3.24. 


Circumferential    Strain    Comparison 
Unstiffened  and  One-Stiffener  Models 


for 
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(a)  Unstif fened  model 
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(b)  Two-stif fener  model 


Figure  3.25 


Longitudinal  Strains  (Unstiffened  and  Two- 
Stiffened)  at  Stiffener  Location  Closest  to 
Charge 
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(a)    Unstiffened  model 
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(b)  Two-stif fener  model 


Figure  3.26. 


Longitudinal  Strains  (Unstiffened  and  Two- 
Stiffened)  at  Stiffener  Location  Far  from 
Charge 
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IV.  CONCLUSIONS 

The  numerical  and  experimental  study  was  performed  to 
investigate  the  nonlinear  dynamic  response  of  cylindrical 
shells  to  end-on  explosions.  Unstiffened  and  ring-stiffened 
cylinders  capped  at  both  ends  were  subjected  to  near-field  and 
far-field  explosions.  Both  finite  element  and  boundary  element 
methods  were  used  for  the  numerical  study.  The  experiment  was 
fulfilled  using  a  6061-T6  aluminum  alloy  cylinder  and  a  sixty 
pound  HBX-i charge  with  28  feet  stand-off  distance. 

The  far-field  explosion  resulted  in  a  nearly  symmetric 
deformed  shape  of  the  cylinder  about  the  center  plane  between 
the  two  end  plates,  while  the  near-field  explosion  caused  very 
unsymmetric  deformation  about  the  center  plane.  Both  far-field 
and  near-field  explosions  induced  localized  deformations  near 
the  end  plates.  The  localized  deformations  were  located  close 
to  both  end  plates  for  the  far-field  explosion  and  a  very 
severe,  localized  deformation  occurred  close  to  the  closest 
end  plate  to  the  explosive  charge  for  the  near-field  explosion 
regardless  of  the  existence  of  ring  stiffeners.  The  accordion 
mode  was  observed  for  the  cylinder  subjected  to  the  far-field 
explosion  but  not  for  the  cylinders  subjected  to  near-field 
explosion . 

A  parametric  study  was  undertaken  to  examine  the  effect  of 
end  plates  on  the  deformation  of  the  cylinder  subject  to  the 
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far-field  explosion.  A  variation  of  the  density  and  stiffness 
of  the  end  plates  caused  a  significant  change  in  the  stress 
wave  propagation  and  the  deformation  of  the  cylinder.  However, 
the  change  of  stiffness  had  more  significant  effects  than  the 
change  of  density. 

The  circumferential  strains  under  the  far-field  explosion 
had  large  initial  compressive  peaks  and  returned  to  the  steady 
state  values  quickly.  The  circumferential  strains  were  larger 
near  the  end  plates  than  around  the  center  due  to  the  local 
deformation.  The  localized  deformations  were  caused  by  the 
compressive  shock  pressure  applied  on  the  end  plates.  The 
longitudinal  strains  under  the  far-field  explosion  were 
compressive  at  the  shell  near  the  end  plate  closer  to  the 
charge  but  tensile  at  the  opposite  location. 

The  failure  of  pressure  and  several  strain  gages  prevented 
a  quantitative  comparison  between  the  numerical  and 
experimental  results.  A  qualitative  comparison  was  possible 
between  the  two  solutions  only  at  a  few  locations.  The 
comparison  was  better  near  the  remote  end  plate  than  the  near 
end  plate  to  the  charge.  A  good  explanation  for  this 
discrepancy  was  not  possible  due  to  the  lack  of  information  in 
the  experiment.  However,  this  was  the  first  test  among  a 
series  of  experiments  to  be  performed.  More  detailed 
information  from  the  next  tests  may  provide  a  better 
understanding  of  the  discrepancy. 
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The  effect  of  ring  stiffeners  was  larger  on  the 
circumferential  strain  than  on  the  longitudinal  strain.  The 
ring  stiffener  located  close  to  the  severe  local  deformation 
had  a  more  pronounced  effect  on  the  strains.  The  ring 
stiffeners  in  general  reduced  the  circumferential  strain 
compared  to  the  unstiffened  case  and  stiffeners  also  gave  more 
recovery  of  the  strain.  However,  when  the  stiffeners  had  a 
large  plastic  deformation,  they  allowed  less  relief  of  the 
circumferential  strain. 
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APPENDIX:  UNIAXIAL  TENSION  TEST  DATA  FOR  6061-T6  ALUMINUM 

A  uniaxial  tension  test  was  performed  on  a  test  specimen 
of  6061-T6  Aluminum  to  find  its  material  properties.  The 
specimen  was  cut  to  ASTM  E8-69  specifications  for  rectangular 
tension  test  specimens. 

A  Material  Test  System  (MTS)  model  810  was  used  for  the 
test  and  the  gage  length  was  set  at  two  inches  for  the  tensile 
test.  The  tensile  test  was  performed  at  a  rate  of  400  sec/inch 
and  a  lad  versus  displacement  curve  was  plotted  as  shown  in 
Figure  A.l.  The  MTS  model  did  not  have  data  recording 
capabilities  and  therefore,  some  points  were  taken  from  the 
graph  after  the  test  to  plot  the  appropriate  stress-strain 
curve.  Some  sample  data  points  are  listed  in  Table  A.l. 

From  the  stress-strain  curve  in  Figure  A.  2,  the  yield 
stress  was  found  to  be  close  to  43  ksi  using  a  0.2  percent 
offset  and  the  modulus  of  elasticity  was  10800  ksi. 
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Figure  A. 2.  Stress-Strain  Curve  for  6061-T6  Aluminum 
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TABLE 

A.l 

LOAD  (LBF)      STRESS 

(psi) 

STRAIN 

0.0 

0.0 

0.0 

500 

10526 

0.0010 

625 

13158 

0.00125 

1000 

21053 

0.0020 

1250 

26316 

0.0025 

1500 

31579 

0.003125 

1625 

34211 

0.003375 

1750 

36842 

0.003625 

1875 

39474 

0.003875 

2000 

42105 

0.004625 

2075 

43684 

0.006625 

2100 

44210 

0.008125 
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